Quantum noise limited and entanglement-assisted magnetometry 
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We study experimentally the fundamental limits of sensitivity of an atomic radio-frequency mag- 
netometer. First we apply an optimal sequence of state preparation, evolution, and the back-action 
evading measurement to achieve a nearly projection noise limited sensitivity. We furthermore ex- 
perimentally demonstrate that Einstein-Podolsky-Rosen (EPR) entanglement of atoms generated 
by a measurement enhances the sensitivity to pulsed magnetic fields. We demonstrate this quantum 
limited sensing in a magnetometer utilizing a truly macroscopic ensemble of 1.5 • 10 12 atoms which 
allows us to achieve sub-femtoTesla/\/Hz sensitivity. 



Ultra-sensitive atomic magnetometry is based on the 
measurement of the polarization rotation of light trans- 
mitted through an ensemble of atoms placed in the mag- 
netic field [1] . For Na atoms in a state with the magnetic 
quantum number rap = F along a quantization axis x 
the collective magnetic moment (spin) of the ensemble 
has the length J x = FN a- A magnetic field along the 
y axis causes a rotation of J in the x — z plane. Po- 
larization of light propagating along z will be rotated 
proportional to J z due to the Faraday effect. From a 
quantum mechanical point of view, this measurement is 
limited by quantum fluctuations (shot noise) of light, the 
projection noise (PN) of atoms, and the quantum back- 
action noise of light onto atoms. PN originates from 
the Heisenberg uncertainty relation SJ Z ■ 5J y > J x /2, 
and corresponds to the minimal transverse spin vari- 
ances SJzy = J x /2 — FNa/2 for uncorrelated atoms 
in a coherent spin state [2]. Quantum entanglement 
leads to the reduction of the atomic noise below PN and 
hence is capable of enhancing the sensitivity of metrol- 
ogy and sensing as discussed theoretically in [2j[9]. In 
[TOT ITT] entanglement of a few ions have been used in 
spectroscopy. Recently proof-of-principle measurements 
with larger atomic ensembles, which go beyond the PN 
limit have been implemented in interferometry with 10 3 
atoms [T2j , in Ramsey spectroscopy [Ml [14] with up to 
10 5 atoms, and in Faraday spectroscopy with 10 6 spin 
polarized cold atoms [15] . 

In this Letter we demonstrate PN limited and 
entanglement-assisted measurement of a radio-frequency 
(RF) magnetic field by an atomic caesium vapour magne- 
tometer. In the magnetometer J precesses at the Larmor 
frequency Q/2tt = 322kHz around a dc field B = 0.92G 
applied along the x axis and an RF field with the fre- 
quency fl is applied in the y — z plane (Fig. [iji). The 
magnetometer (Fig. [2^,) detects an RF pulse with a con- 
stant amplitude -Brp and duration r (Fourier limited full 



width half maximum bandwidth S = 0.88r 



)• 



The mean value of the projection of the atomic spin on 
the y — z plane in the rotating frame after the RF pulse 
is TB RF J X T 2 [1 - exp(-r/T 2 )]/2. Here T 2 is the spin de- 
coherence time during the RF pulse and V = fl/B = 



2.2 • 10 10 rad/(sec-Tesla) for caesium. Equating the mean 
value to the PN uncertainty we get for the minimal de- 
tectable field under the PN limited measurement 
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The PN-limited sensitivity to the magnetic field is then 
-BminV^j an( i equals the standard deviation of the mea- 
sured magnetic field which can be achieved by using re- 
peated measurements with a total duration of 1 second. 
The best sensitivity to the Brp with a given S is achieved 
with the narrow atomic bandwidth: 5 = t _1 > T 2 ~ l ■ 
A long T 2 also helps to take advantage of the entangle- 
ment of atoms. Entangled states are fragile and have a 
shorter lifetime T 2 e < T 2 . As we demonstrate here, un- 
der the condition (5 _1 < T 2 e < T 2 , that is for broadband 
RF fields, entanglement improves the sensitivity. Similar 
conclusions have been reached theoretically for atomic 
clocks in [3l [5] and for dc magnetometry in [4] . 

In our experiment a long magnetometer coherence time 
T 2 > 30msec is achieved by using paraffin coated caesium 
cells at around room temperature [16] , and by the time 
resolved quantum spin state preparation (optical pump- 
ing), evolution, and measurement. In this way T 2 is not 
reduced by the optical pumping and/or measurement- 
induced decoherence during the time when magnetic field 
is applied. 

PN limited sensitivity requires, besides elimination of 
the technical noise, the reduction of the back action noise 
of the measurement and of the shot noise of the probe 
light. 

The back action noise comes from the Stark shift im- 
posed by quantum fluctuations of light polarization on 
J y ab in the laboratory frame when J l z ab is measured. 
When atoms are exposed to an oscillating magnetic field 
J l z ab and J l y h experience Larmor precession, and hence 
both of them accumulate the back action noise [T7] [TB] . 
As proposed in [T7] the effect of this noise can be canceled 
if two atomic ensembles with orientations J x and — J x are 
used. Fig. [2^l presents the sketch of the magnetometer 
layout which is used for the back action evading measure- 
ment of the Brf , and Fig. [3] shows the level schemes for 
the two atomic ensembles. 
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FIG. 1: (a) Radio-frequency magnetometer. The atomic spin 
J precesses in crossed dc and RF magnetic fields (the blue 
dashed spiral). The precessing J imposes an oscillating po- 
larization rotation on the probe light, (b) Results of a series 
of measurements of S2s,c without Bef (g re y points) and af- 
ter a 15 msec RF pulse with Brf = 36 fT (red points). The 
solid and the dashed lines show the standard deviation of the 
experimental points, and of the the projection noise (PN) lim- 
ited measurement, respectively. The spin displacement in the 
rotating frame is indicated by the dashed blue arrow. (Inset) 
Power spectrum of the photocurrent (arb. units, centered at 
57) The narrow peak of the red/black spectrum is the atomic 
spin signal/projection noise. The broad part of the spectrum 
is the shot noise of light with the width set by the detection 
bandwidth. 



To read out the atomic spin, which carries information 
about the RF magnetic field, we utilize polarized light. 
The quantum state of light is described by the Stokes 
operators S±, S2 and S3. The relevant observables are 
the cos/sin f2 Fourier components of S2 integrated over 
the pulse duration T with a suitable exponential mode- 
function §2 C <x Jq S2(t) cos(ftt)e ±7tot *. S\ s is similarly 
defined by replacing cosine with sine. The operators S^c 
and S*2 S are experimentally measured by homodyning and 
lock-in detection (Fig. [2]}) . 

The theory of interaction between the probe light and 
a vapor of spin polarized alkali atoms has been developed 
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FIG. 2: (a) Schematic of the experimental setup. Cubic 22mm 
paraffin coated cells filled with caesium are placed in the dc 
field B in two separate magnetic shields. The atoms are opti- 
cally pumped so that the directions of the collective spins in 
the two cells (black arrows) are opposite. A pulse of Brf at 
the frequency fi is applied orthogonally to the B field. The 
polarization rotation of the top-hat shaped probe beam pulse 
(diameter 21mm) is detected by two detectors (HWP - half 
wave plate). The lock-in amplifier measures the cos fit and 
sin Qt components of the photocurrent. (b) Pulse sequence for 
projection noise limited magnetometry. The temporal mode 
function for the probe is shown with a dashed black curve, (c) 
Pulse sequence and temporal modes for entanglement-assisted 
magnetometry. 












a h i / 


1 \a. 




\ 1 


r-4 


b 2 



FIG. 3: Level scheme and optical fields used in the magne- 
tometer. The input probe light polarized orthogonally to the 
magnetic field (solid arrows) is blue detuned by A = 850MHz 
from the F = 4 — > F = 5 hyperfine component of the D2 
line at 852nm. Forward scattering of photons into modes a± 
leads to polarization rotation and swapping interaction de- 
scribed in the text. The two diagrams correspond to the two 
oppositely polarized cells using which leads to cancelation of 
the backaction of light onto atoms, as explained in the text. 



in [121 [20] ■ After the probe has interacted with the two 
atomic ensembles, the output Stokes operator (normal- 
ized for a coherent state as Var (S) = $/2, with $ being 
the photon number per pulse) is given by [20] 
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with the equation for SJ" ' obtained by the replacement 
J" 1 — > Jy l . The output Stokes operator is defined with 
an exponentially falling mode function e _7tot , and the in- 
put Stokes operator is defined with a rising mode e +7tot * 
[2"U] , where 7 to t is the decay rate of the transverse atomic 
spin components, k — — cxp(—2j swap T) is the 
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light-atom coupling constant, 7 swap is proportional to in- 
tensity of light and density of atoms, and £ = ^/Tia^Jai, 
where ai and a\ are the tensor and vector polarizabilities. 
For our probe detuning A = 850MHz, £ 2 = 1/6.3. In the 
two-cell setup 52c and Bi s contain information about the 
commuting rotating frame operators J*™ = </*™ + J*5 an d 
J % z \ = Jli + J l Z 2 simultaneously, which are displaced by 
Brf- 

The atomic spin operator J z+ after the interaction is 
given by 



jout 
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In case of 7 SW apT *C 1, corresponding to either rather 
large probe detuning A, or to a relatively small pho- 
ton number <&, as in PHUT], Eq. Q and @ reduce to 
the quantum nondemolition (QND) measurement, where 
J y + and J z + are conserved. Here we implement a strong 
measurement limit, 7 SW apT ~ 1 in which the light and 
atoms nearly swap their quantum states. Indeed under 
this condition first terms in Eq. ^) and ^ are strongly 
suppressed. The exponential suppression of the probe 
shot noise contribution to the signal with time and with 
the optical depth d in Eq. ^) allows to approach the PN 
limited sensitivity faster than in the QND measurement 
case where this contribution would stay constant. 

The PN limit can be approached if the total rate of the 
decay of the atomic coherence 7tot is dominated by the 
coherent swap rate 7 swa p. In the experiment we obtained 
7swap = 0.43ms _1 and jtot = 0.50ms _1 with a 15mW 
probe detuned by 850MHz and N A = 2 • 7.2(7) • 10 11 
corresponding to the effective resonant optical depth d — 
75. The residual decoherence rate due to the spontaneous 
emission rate 771, collisions, and magnetic dephasing is 
therefore 0.07ms _1 . Note that 7 swap ~ £ 2 d • 771 3> T]\. 
The probe losses are dominated by reflection on cell walls 
and on detection optics, since its absorption in the gas is 
negligible. 

Raw results for a series of measurements of Sfc s °b~ 
tained at 32°C corresponding to N A = 2 • 7.2(7) • 10 11 
are presented in Fig. [TJd. N a is monitored via J x = 4N A 
measured by the Faraday rotation angle 19° of an aux- 
iliary probe beam sent along the x axis, and from the 
degree of spin polarization > 98%, as determined from 
the magneto-optical resonance [22] • The experimental 
sensitivity is B RF ■ y/r/SNR = 3.6(4) • 10" 16 Tesla/\/Hz, 
where the signal to noise ratio SNR= 12.3 is found from 
the data in Fig. [TJd as the ratio of the mean to the 
standard deviation of the data (red points) obtained for 
Brf = 36(3)fT applied during t = 15msec via a cal- 
ibrated RF coil. The optimal temporal mode function 
for this measurement has j p t = 2j t ot = 1.0ms . This 
sensitivity is ~ 30% above the theoretical PN limited 
sensitivity 2.7(5) • 10- 16 Tesla/VHz found from (Eq. Q) 
using the measured value of T2 = 32ms. Using this sensi- 
tivity we can calibrate the values in Fig. lb in PN units. 




FIG. 4: Entanglement assisted magnetometry results, (a) 
Magnetometer read out noise in units of shot noise of light. 
The dashed line at 0.5 is the contribution of the probe light 
noise(LN). The average noise level of the initial spin state 
(squares) is 1.10(8) in projection noise (PN) units. The read 
out noise corresponding to the PN level is shown as the solid 
horizontal line. Grey bands have the width equal to twice 
the standard deviation. The magnetometer noise for the en- 
tangled state vs the delay between the two probe pulses is 
shown as circles. The fit yields the entanglement lifetime 
T2E = 4msec. (b) Sensitivity of the entanglement assisted 
magnetometer (circles) and the sensitivity in the absence 
of entanglement (squares) vs the RF pulse bandwidth (in- 
verse pulse duration). The Brf = 36 fT pulse is applied to 
Na = 3.6(3) • 10 11 atoms. The duration of the first/second 
probe pulse is 2/3msec, the duration of the pumping cycle 
6msec. (c) Signal/noise ratio times the RF bandwidth for the 
magnetometer in the presence (circles) and in the absence 
(squares) of the entanglement. 



The difference between the two sensitivities is due to the 
residual contribution of the shot noise of the probe (first 
term in Eq. pi), the decay of the atomic state during 
the RF pulse and the classical fluctuations of the atomic 
spins. 

The experimental sensitivity B RF -VT/SNR = 4.2(8)- 
10~ 16 Tesla/vTIz calculated using the full measurement 
cycle time T including the duration of optical pump- 
ing (6ms), probing (1.5ms) and r = 22ms (Fourier lim- 
ited bandwidth ~ 40Hz) approaches the best to-date 
atomic RF magnetometry sensitivity [2 3) obtained with 
10 4 times more atoms. This is to be expected since PN 
limited magnetometry yields the best possible sensitivity 
per atom achievable without entanglement. 

We now turn to the entanglement-assisted magnetom- 
etry. As first demonstrated in |17j the measurement 
of the Stokes operators Sq^\ can generate the state 
of atoms in the two cells which fulfils the Einstein — 
Podolsky — Rosen (EPR) entanglement condition for two 
atomic ensembles with macroscopic spins J x and — J x : 



J EPR 



Var i 



Tin 
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J z2 



/2J X < 1. 



(jyi + Jyi) + Var ( u 
This inequality means that the atomic spin noise which 
enters in Eq. ^ is suppressed below the PN level corre- 
sponding to the coherent spin state. Entanglement can 
be visualized as correlations between the grey and red 
points shown in Fig. [T]d. The degree of entanglement 
and the sensitivity are optimized by choosing suitable ris- 
ing/falling probe modes e ±7 ' for the first/second probe 
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modes shown in (Fig. |2p). In order to reduce the contri- 
bution of the technical (classical) fluctuations of the spins 
we perform the entanglement-assisted measurements at 
room temperature with Na = 3.6 • 10 11 per cell. 

In order to demonstrate entanglement of atoms, we 
need to calibrate the PN level. Above we have calcu- 
lated the atomic noise in PN units with 20% uncertainty 
using (Eq. Q). However this uncertainty coming from 
the uncertainties of Na and Brf is too high and we 
therefore apply the following calibration method. The 
collective spin operator = J y i — J y2 after the in- 
teraction with a probe pulse is given by [20]: J°"* = 

{Jyl ~ jfflV 1 -^* 2 + K^^w^Sic- Using an electro- 
optical modulator we create a certain average value of 
S>1™ and then calibrate it in units of shot noise -\/$/2 by 
homodyne detection. We determine k 2 for a certain light 
power P « 8mW by first sending a 1msec pulse with 
(S^™) = y / $/2, which is mapped on the atoms. We then 
flip the sign of J y2 and J z2 by momentarily increasing the 
dc magnetic field and send a second 1msec light pulse. 
Now (J v -) is mapped on (S^), following Eq. and 
the measurement on the second pulse yields = n 2 . 

Using k 2 and the detection efficiency r\ = 0.8 we can cal- 
culate the atomic noise in projection noise units from the 
measured noise Var yS^sj usm g Eq. With a known 
value of k, we can now calibrate the atomic displace- 
ment {(Jy), (J z )) caused by a particular Brp in units of 
PN as follows. The mean value of the homodyne signal 
is (k ut ) = n^<S>/2FN A {{J yl + J y2 ) + i{J zl - J z2 )). 
From this expression the atomic displacements in units 
of PN \JFNa can be found using the values of k, rj and 
the shot noise $/2. Once we know the atomic displace- 
ment in units of PN for the RF pulse with a certain du- 
ration and amplitude we can utilize this to find n and 
hence the PN level for probe pulses with any light power 
(most of the measurements were done with a probe power 
P w 14mW) from the mean value of the homodyne sig- 
nal corresponding to the RF pulse with the same duration 
and amplitude. 

Entanglement is generated in our magnetometer by the 
entangling light pulse applied before the RF pulse in the 
two-pulse time sequence (Fig. [2J:). Fig. [4^, shows the 

magnetometer noise Var (^S 2 wt ^j including the light noise 
and the atomic noise contributions in units of shot noise 
of the 15mW probe, as a function of the RF pulse du- 
ration for entangled atoms and for atoms in the initial 
state. For these measurements k 2 — 3.1 which is con- 
sistent with Fig. [4ji since, according to Eq. for PN 
limited measurement k 2 is equal to the total noise 3.6 less 
the light noise 0.5. Knowing k 2 we find the initial noise 
level in PN units to be 1.10(8). The best noise reduction 
below PN of about — 1.5dB (-30%) is obtained for short 
RF pulses corresponding to the RF bandwidth 5rf > 1 
kHz. Fig. [4Jd illustrates the improvement in the sensitiv- 



ity Brf • VT/SNR with entangled atoms compared to 
the sensitivity obtained with atoms in the initial state. 
Fig. [4}: shows that entanglement improves the product of 
the S /N ratio times the bandwidth Srf by a constant fac- 
tor for RF pulses with Srf much greater than the inverse 
entanglement lifetime (itT2e)~ 1 — 70Hz. 

To the best of our knowledge, our results present en- 
tanglement assisted metrology with the highest to-date 
number of atoms. This results in the absolute sensitiv- 
ity at the sub-femtoTesla/vTIz level which is comparable 
to the sensitivity of the state-of-the-art atomic magne- 
tometer |23| operating with orders of magnitude more 
atoms. The two-cell setup can also serve as the magnetic 
gradient-meter if the direction of the probe in the sec- 
ond cell is flipped or if the RF pulse is applied to one cell 
only. Increasing the size of the cells to a 5 cm cube should 
yield the sensitivity of the order of 5 • 10 _17 T/\/Hz which 
approaches the sensitivity of the best superconducting 
SQUID magnetometer [Mj. The degree of entanglement 
Sepr can, in principle, reach the ratio of the tensor to 
vector polarizabilities £ 2 which is 0.16 or — 8dB for our 
experiment, and can be even higher for a farther detuned 
probe. This limit has not been achieved in the present 
experiment due to various decoherence effects, including 
spontaneous emission and collisions. Increasing the size 
of the cells may help to reduce some of those effects since 
a larger optical depth will then be achieved for a given 
density of atoms. 
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